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Using the generalized gradient approximation augmented with maximally localized Wannier functions 
analysis, we present the formation of cuprate-like electronic structures in Ag n F 2 -related superlattices 
resulted from the confinement together with structural chemical modification. The out-of-plane electronic 
reconstruction leads to electron doping of AgF 2 plane and gradually destablizes the antiferromagnetic state. 
Eventually a stable nonmagnetic metallic state emerges by applying in-plane tensile strain, in which the 
shape of effective Fermi surface of AgF 2 plane exhibits the key feature of high-temperature cuprate 
superconductor. 

The discovery of high-temperature superconductivity in cuprates 1 initiated the quest for exploring super- 
conductivity in other transition-metal compounds. Ag 2+ is isoelectronic with Cu 2+ (d 9 configuration). 
F~ and O 2- are also isoelectronic ions, closed-shell species, moreover both F" and O 2 " are weak- field 
ligands. Given the existence of the superconducting cuprates, one might be naturally interested in searching the 
superconductivity in Ag 2+ -F~ solids. Theoretical study 2 " 4 of Grochala and Hoffmann has also suggested that 
properly hole- or electron-doped Ag n fluorides might be good superconductors, due to similarity in structure and 
properties between the Ag n fluorides and the cuprate superconductors. Particularly, Jaron and Grochala have 
predicted that the high pressure a form of AgF 2 compound (>15 GPa) is a layered material with antiferromag- 
netic (AFM) ordering 4 , as seen for the parent compounds of high-temperature cuprate superconductors 5 . 

It is well known that [Cu0 2 ]oo plane with tetragonal tetra coordination of Cu (weak apical Cu-O bonds), is an 
essential structural element for superconductivity in cuprates. The basic band structure of the doped cuprates is a 
single 2D C\x-d x i _ y i — like band deviating from half filled. In this situation, AFM fluctuations prevail in the 
undoped parent compounds and are often believed to mediate the superconductivity. The Fermi surface (FS) 
from this d X 2- y2 band has been observed in many overdoped cuprates and agrees with the density functional 
theory based band structure calculations 5,6 . 

As Grochala and Hoffmann pointed out in their review paper 2 , analogous [AgF 2 ]o° plane with tetracoordinate 
Ag has not been reported experimentally. Noteworthy recent development of heterostructure interface techno- 
logy, superlattices (SLs) containing Ag n F 2 square lattices can be prepared by using appropriate synthetic tech- 
niques to incorporate alternating layers of different transition metal compounds 7 " 14 , even technically a single 
atomic layer 11 . Here, interface can be used to modulate electronic structure for manipulating physical properties 
and generating novel phases which are not present in the bulk constituents. Whereby the quest for Ag 11 super- 
conductors could be achieved by the aforementioned novel paradigm on designing and fabricating two dimen- 
sional materials 15 " 18 . In our paper, our research focuses on artificial superlattice materials design and their 
electronic properties and doping effects, different from research on real bulk Ag 11 fluorides materials 2 " 4 . 

To identify possible superconductivity in 2D Ag 2+ -F" square sheet, we investigate electronic structures, 
magnetic states, model hamiltonian parameters and effective FSs for three proposed superlattices: SrTi0 3 / 
AgF 2 (STO/AgF 2 ), BaTi0 3 /AgF 2 (BTO/AgF 2 ) and SrTi0 3 /CsAgF 3 (STO/CsAgF 3 ), as illustrated in the top panels 
of Fig. 1, and compare these with corresponding properties of the cuprate superconductors. We find cuprate-like 
band structures and strong AFM fluctuations in all these SLs. More importantly, large cation size increases 
cation-anion polarization strength and corresponding apical Ag-O/F distance, and makes oxygen band edge 
shift above the Fermi level to exchange charge with Ag-d x 2 _ y i band, leading to out-of-plane electronic recon- 
struction. Here, similar to charge transfer in cuprates and recently studied La 2 Cu0 4 -related heterostructure 18 , the 
Ti0 2 plane actually serves as a charge reservoir block and transfer electrons to the AgF 2 plane. As a result, AFM 
state presents an unstable trend. Finally, the applied in-plane tensile strain drives a novel phase transition from 
the AFM metallic state to a stable nonmagnetic (NM) metallic state in STO/CsAgF 3 SL. Model hamiltonian 
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Figure 1 | Schematic geometrical structures and GGA bandstructures of bulk LCO, bulk HBCO, STO/AgF 2 , BTO/AgF 2 and STO/CsAgF 3 SLs from left 
to right. The Fermi level e F is set at zero. Dark cyan and orange fatbands represent contribution of d x 2_ y i and d 3z 2_i orbitals respectively. 



parameters and FS character of STO/CsAgF 3 are extracted and com- 
pared to La 2 Cu0 4 (LCO) and HgBa 2 Cu0 4 (HBCO), indicating that 
STO/CsAgF 3 is a promising candidate for Ag n superconductivity. 

Results 

For the in-plane lattice constant a of SLs, we first took that of STO 
(3.905 A), often used as the substrate. The lattice constant c and 
atomic z coordinates were fully relaxed. The main effect of the relaxa- 
tion, is to make the negatively charged O/F and positively charged 
cations displaced relative to each other in SrO, BaO and CsF atomic 
layers, and thereby polarize the cation and anion planes. Cation size 
affects significantly polarization strength and corresponding apical 
Ag-O(F) bond length. AgF 2 layer acts as the mirror plane of whole 
unit cell. In STO/AgF 2 and BTO/AgF 2 SLs, oxygen atoms move 
symmetrically towards and against AgF 2 plane by 0.077 A and 
0.064 A respectively. As a result, apical Ag-O bond length in the 
BTO/AgF 2 is slightly larger than that in the STO/AgF 2 by 0.143 A, 
due to the larger size of the Ba 2+ cation. The largest cation-anion 
polarization occurs in CsF plane in the STO/CsAgF 3 , and is 0.623 A 
against AgF 2 plane. This polarization distortion produces a local 
ionic dipole moment, and together with in-plane strain it also leads 
to a larger increment in the apical Ag-F distance (d a £ g lc f F = 3.405 A) 

compared to apical Ag-O distance (d^f 0 = 2.486, 2.629 A) in 
other two SLs. 

Local ionic dipole moment perturbs electrostatic potential and 
changes band positions around the Fermi Level. An evolution of 
Ag-e g states with structural chemical modification can be clearly 



observed in band structures. Spin-polarized GGA calculations give 
paramagnetic ground state for all these superlattices. Fig. 1 shows 
energy bands of STO/AgF 2 , BTO/AgF 2 and STO/CsAgF 3 SLs in a 
13 eV region around the Fermi level s F = 0 and along the symmetry- 
lines r-X-M-r = (0,0,0) - (- ,0,o) - (- , - ,0] - (0,0,0) . The 

energy bands of bulk LCO and HBCO are also plotted in Fig. 1 for 
comparison. For three SLs, electronic properties around e F are still 
mainly controlled by Ag-e g bands, which are above the filled OIV-2p 
and Ag-t 2g bands, and below the empty Ti-3d bands. We plot d x i _ y i 
(dark cyan) and d^ z i _ i (orange) fatbands around s F to disclose their 
orbital contribution. Ag-e g antibonding bands have similar band 
width in superlattice configuration STO/AgF 2 and BTO/AgF 2 , and 
resemble that of LCO. Note that oxygen 2p bands become closer to 
Fermi Level in BTO/AgF 2 due to atomic polarization distortion in z 
direction. However, in STO/CsAgF 3 case, the antibonding band 
between Ag-d 3z 2 _ x and F-p states disappears due to the weak mixing 
of Ag-3d and F-p states in z direction. As a result, e g bands from — 4 to 
2 eV looks more like that of HBCO with a larger apical Cu-O dis- 
tance of 2.784 A. Most importantly, oxygen 2p band edge of Ti0 2 
plane upshifts eventually above the Fermi level and exchanges charge 
with Ag-d x 2_ y 2 band, as seen in layer-projected density of states in 
the left panel of Fig. 2. 

In order to investigate the microscopic orbital physics of polariza- 
tion-induced electron doping, we plot the typical partial charge den- 
sity of the unoccupied bands of STO/CsAgF 3 from Fermi level to 
0.25 eV in the right panel of Fig. 2. Obviously, electron doping 
originates from charge transfer between 0-p x ,p y orbitals at Ti0 2 
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Figure 2 | Layer-projected density of states (left panel) of STO/CsAgF 3 SL, and corresponding partial charge density isosurfaces (right panel) for the 
unoccupied bands from Fermi level to 0.25 eV. 



interface and Ag-d x 2 _ y i orbital. Here, Ag-<i x 2 _ y i has considerable 
covalent hybridization with in-plane fluorine atoms' p x , p y orbitals, 
similar to Cu 2+ -0 2 ~ bonds. Interestingly, 0-p x ,p y state around the 
Fermi level is mainly located at Ti0 2 interface and far away from the 
AgF 2 plane, as illustrated in Fig. 2, which favors the realization of 
possible superconductivity in 2D Ag n F 2 plane. 

Next, we discuss the stability of magnetic states in three super- 
lattices under GGA + U d scheme. AFM band structures with Ag-e g 
orbitals character (see Fig. 1 of supplementary materials) indicate 



that STO/AgF 2 SL presents an AFM insulating ground state with a 
energy gap of 0.45 eV. While in BTO/AgF 2 electronic correlation 
drives a weak electron doping in Ag-d x 2_ y i state, which is absent 
in GGA bands, leading to an insulator-metal transition. For STO/ 
CsAgF 3 , an AFM metallic ground state is obtained. By analyzing 
layer-projected density of states of STO/CsAgF 3 from GGA + U d 
calculations and corresponding partial charge density isosurfaces for 
the bands around Fermi level (see Fig. 2 of supplementary materials), 
we find that the obtained AFM metallic ground state is aroused by 



Table 1 | The in-plane and apical bond length d™'^' 
atom's magnetic moment of AFM state in /ie/Ag(Cu), for 
and with (bottom subtable) in-plane strain 


and ^Ag(^u\-F(o)} n A, energy differences Enm-Eafm in meV/Ag(Cu), and Ag/Cu 
LCO, f%CO, STO/AgF 2 , BTO/AgF 2 and STO/CsAgF 3 without (top subtable) 






LCO HBCO 


STO 
AgF 2 


BTO 

AgF~ 2 


STO 
CsAgF 3 


jin —plane 
a Ag(Cu)-F(0) 
japical 

a Ag(Cu)-F(0) 

Enm-Eafm 
Moment 




1.894 1.941 

2.429 2.784 

206.510 117.470 
0.542 0.495 


1.953 

2.486 

91.160 
0.513 


1.953 

2.629 

34.525 
0.403 


1.953 

3.405 

4.105 
0.184 










STO 












CsAgF 3 




jin-plane 
a Ag-F 




1.895 




2.003 


2.053 


iapical 
d Ag-F 




3.760 




3.368 


3.289 


Enm-Eafm 
Moment 




7.595 
0.212 




0.335 
0.102 


* 
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*The instability occurs 


to the AFM state. 
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Figure 3 | (a) Parameters of the six-band p-d model for the CuO and AgF(O) octahedral in cuprate superconductors and the proposed AgF 2 -related 
superlattices; (b) Localized Wannier functions of Ag — d x 2_ y 2 and — d 3z i_ 1 orbitals in STO/CsAgF 3 ; (c) Effective Fermi surface centered at T point 
in first Brillouin Zone from d x 2_ f band for bulk LCO, bulk HBCO, STO/AgF 2 , BTO/AgF 2 and STO/CsAgF 3 SLs. 



charge transfer between 0-p x ,p y orbitals in Ti0 2 plane and covalent 
hybrid orbitals of Ag-d x 2 _ y i and ¥-p x ,p y in AgF 2 plane. In all three 
superlattices, FM state falls in between AFM and NM states in 
energy. In Table I, we summarize in-plane and apical bond lengths 



tin —plane 



and d 1 



^apical 



-A g {cu)-F{0) ^ "A g (Cu)-F(oy ener SY difference E NM -E AFM , and 
magnetic moment on Ag/Cu atom in AFM state. With the 
increasing apical Ag-O/F distance, E NM -E AFM value decreases 
gradually from 91.160 meV/Ag for STO/AgF 2 to 34.525 meV/Ag 
for BTO/AgF 2 , similar to the trend for cuprates (e.g. from 
206.51 meV/Cu for LCO to 117.47 meV/Cu for HBCO in 
Table I), and finally to a much smaller value 4.105 meV/Ag in 
STO/CsAgF 3 . The process of AFM state instability is also compa- 
nied by reduction of magnetic moment on Ag/Cu atoms. For 
superlattice structures, electron doping of AgF 2 plane emerges 
with the change of apical Ag-O/F distance, as we discussed above, 
which is an important derivation of AFM state instability, espe- 
cially in STO/CsAgF 3 with much smaller E NM -E AFM value and 
magnetic moment of 0.184 ^/Ag. 

Furthermore, we investigate the effect of in-plane strain on mag- 
netic state, because electronic properties are subject to electron- and 
orbital-lattice couplings in perovskite-like materials. Similar calcu- 
lations are made for STO/CsAgF 3 with three additional in-plane 



lattice constants of 3.790, 4.005 and 4.105 A. We find that compress- 
ive strain can effectively increase E NM -E AFM to 7.595 meV/Ag, but 
tensile strain decreases it to 0.335 meV/Ag for a = 4.005 A. 
Extraordinarily, STO/CsAgF 3 SL goes through phase transition to 
a stable NM metallic ground state under tensile strain a = 4.105 A, 
suggesting that the tuning of in-plane lattice constant can serve as an 
effective tool to modulate magnetic properties and even 
superconductivity. 

We know that effective low- energy hamiltonians containing the 
minimal set of bands are important tools for understanding chem- 
ical trends. Based on the aboved GGA simulations, we extract 
model hamiltonian parameters by MLWFs downfolding technique. 
In this work, we choose to downfold to a 6 -band hamiltonian 
describing the in-plane d x 2_ y i, p x , p y orbitals, and out-of-plane 
d 3z 2_ ly two p z orbitals (see Fig. 3a). In particular, six parameters 
capture the essential physics: the e g crystal field splitting energy 
Ai = 8 X 2 _ y i — g 3z 2 _ l , the in-plane charge-transfer energy 
A 2 = s x 2 _ y 2 — s Px{y) , the direct in-plane and out-of-plane Ag-F 
(Cu-O) hopping t p d and t* d , and the two shortest-ranged O-O 
hoppings tppy and tp p >. The extracted values are tabulated in 
Table II, and corresponding interpolated band structure are shown 
in Fig. 3 of supplementary materials. 
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Table II | Tight-binding parameters of the six-band p-d model, containing the in-plane d x 2_ y z, p Xf p x orbitals and out-of-plane d 3z i_ 1 , p z 
orbitals for LCO, HBCO, STO/AgF 2 , BTO/AgF 2 and STO/CsAgF 3 without (top subtable) and with (bottom subtable) in-plane strain. 



Parameters include e g crystal field splitting energies Ai = 



charge-transfer energies A 2 



, the three nearest- 



neighbor (intra— cell) hoppings t pc j, f d , t pp , and the inter-cell oxygen-oxygen hopping t pp >. The in-plane and apical bond lengths 



iin— plane j japical 

a Ag(Cu)-F(0) anQ a Ag(Cu)- 



jin —plane 



HO) 



.. l pd . , 
in A are also listed 



"Ag(Cu)-F(0) 
japical 

a Ag(Cu)-F(0) 

Ai 
A 2 
t p d 

t z pd 



LCO 

1.894 

2.429 

0.005 
2.305 
1.395 
0.458 
0.656 
0.191 



HBCO 

1.941 

2.784 

0.115 
1.476 
1.249 
0.391 
0.620 
0.225 



STO 

1.953 

2.486 

0.268 
3.703 
1.766 
0.972 
0.391 
0.089 



BTO 

Wi 

1.953 

2.629 

0.324 
3.729 
1.759 
0.830 
0.398 
0.091 



STO 
CsAgF 3 
1.953 
3.405 

0.810 
3.504 
1.755 
0.223 
0.432 
0.118 



Ai 
A 2 



STO 
CsAgF 3 


1.895 


2.003 


2.053 


3.760 


3.368 


3.289 


1.098 


0.654 


0.520 


3.968 


3.122 


2.785 


1.979 


1.583 


1.431 


0.141 


0.227 


0.245 


0.479 


0.397 


0.364 


0.141 


0.104 


0.092 



The hopping integrals t pd and t pp of LCO and HBCO are in good 
agreement with the 3 -band model results by Weber et al 30 . While A 2 
and tpp' are further corrected in our model by including three addi- 
tional out-of-plane orbitals. From Table II, we find that cuprates and 
Ag n F 2 -related SLs share some common features. The larger d^f F ^ 

or d U Q™ l 0 value leads to the increasing A x and decreasing out-of- 
plane hopping f d respectively. And A 2 , in-plane hoppings t pd and t pp 

increase with the decreasing d l ^~J l F ane (see bottom subtable) or 

dcu-o™ respectively. However, in-plane hopping t pp > is an excep- 
tion, and is affected considerably by out-of-plane Cu-O distance. 
From LCO to HBCO, the weakened electrostatic repulsion enhanced 
the hopping t pp > by 0.034 eV when the negatively- charged apical 

oxygen moves against the CuO plane, although d™']* 1 ™ 6 increases 
from 1.894 to 1.941 A at the same time. The similar feature can be 
seen more clearly in three SLs with the fixed d^~J F ne of 1.953 A. t pp > 
changes from 0.089 eV to 0.091 eV, and finally 0.118 eV with the 
increasing apical Ag-O/F distance. Compared to cuprates, generally 
STO/CsAgF 3 has relatively larger A l5 A 2 , and in-plane Ag-F hopping 
tpd, while hopping f pd , t pp and t p p> are smaller. After applying in-plane 

tensile strain d l ^~^ F ane = 2.003, 2.053 A (see bottom subtable), one 
can find that parameters evolve towards those in cuprates, except for 
tpp and tpp'. 

In Fig. 3b, we plot the localized Wannier functions of Ag-d x 2 _ y i 
and — d 3z 2 _ l in STO/CsAgF 3 SL. Both look like those of cuprates, 
and have a strong p-d covalent hybridization characteristic. 
Ag-d 3z 2_ 1 is more localized due to the big apical Ag-F distance 
(>3 A). The FSs centered at T point for LCO and HBCO are shown 
in the first row of Fig. 3c. Compared to LCO (transition temperature 
T c = 40 K), the FS of HBCO (T c = 90 K) has the typical shape of 
high-T c cuprates superconductor with constant-energy surface 
obviously bulging toward T point. The FS shape of STO/AgF 2 and 
BTO/AgF 2 (see the second row of Fig. 3c) is far away from that of 
HBCO or LCO. However, for STO/CsAgF 3 (the third row of Fig. 3c) 



with polarized electron-doping in AgF 2 plane, effective FSs from 
Ag-d x 2_ y 2 band presents the considerable similarity to that of 
HBCO. 

Discussion 

We analyze the cuprate-like electronic structures and strong AFM 
fluctuations in the proposed superlattice with 2D Ag n F 2 square sheet. 
Atomic polarization induces out-of-plane electronic reconstruction 
occurring between 0-p x ,p y orbitals in Ti0 2 plane and covalent 
hybrid orbitals of Ag-d x 2 _ y 2 and ¥-p x ,p y in AgF 2 plane, which is an 
important origin of AFM state instability. A stable NM metallic 
ground state emerges in STO/CsAgF 3 SL subjected to in-plane 
tensile strain, meanwhile corresponding Wannier functions of 
Ag-d x 2_ y 2, Ag-d 3z 2_ u and FS shape present considerable similarity 
to those in cuprates with the high-T c . Therefore, d 9 Ag n F 2 -related 
superlattices are promising because their physics contains the main 
ingredients of high-temperature superconductivity. 

Method 

We carried out the numerical calculations using the Vienna ab initio Simulation 
Package (VASP) 19-22 within the framework of the generalized gradient approximation 
(GGA) (Perdew-Burke-Ernzerhof exchange correlation functional) 23 , and recently 
developed maximally localized Wannier functions (MLWFs) downfolding tech- 
nique 24 " 26 . The ion-electron interaction was modeled by the projector augmented 
wave (PAW) method 27,28 with a uniform energy cutoff of 500 eV. Spacing between k 
points was 0.02 A" 1 . The structures of the SLs were optimized by employing the 
conjugate gradient technique, and in the final geometry, no force on the atoms 
exceeded 0.01 eV/A. For magnetic states calculations, we used U d = 7.5 eV and J d = 
0.98 eV for Ag-d and Cu-d states 29 . 
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